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1. Introduction 

Uniaxially modulated structures are observed in very different classes of magnetic and 
ferroelectric substances. In many cases they exhibit rather complex phase diagrams 
with large varieties of phases. Phase transitions from a high temperature paramagnetic 
or paraelectric phase (paraphase) to commensurately and incommensurately modulated 
phases occur, as external control parameters like temperature and elastic stresses are 
varied. Microscopic models are successfully used for the description of these modulated 
systems. They were reviewed e.g. in |T|. An interesting example is the p-state chiral 
clock model [H, whose Hamiltonian is 



H = -JqYj H cos 
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a labels the layers perpendicular to the direction of the modulation (chiral direction) 
and i, j the crystal units in these layers, (ij) runs over neighbouring pairs in the layers. 
The integer variables nj^„ describe the state of the unit {i,a). They assume one of the 
values from to p — 1. Below they are called spins. The two terms in equation (|l|) 
describe couplings (Jo > 0, J > 0) between nearest neighbours in the same and in 
adjacent layers, respectively. 

In the ground state every layer is ferromagnetically ordered. Depending on the value of 
A, various ordering patterns of the different layers are realized. For < A < ^ nearest 
neighbours in the chiral direction couple ferromagnetically (ferromagnetic bond), thus 
leading to a ferromagnetic ground state where all spins are equal. For | < A < 1 the 
spin increases by one for successive layers (chiral bond), thus yielding the right-handed 
chiral pattern 

... 012 ... (p - 1) 1 ... 

A = I is a multiphase point at which an infinity of different phases are degenerate since 
ferromagnetic and chiral bonds have the same energy. 

Whereas the three-state model (p = 3) f^, |^ has been very thoroughly investigated, 
only few results are known for the general case p > 4. There are derivations of the 
low temperature phase diagram of the general p-state model by an expansion of the 
free energy in the vicinity of the multiphase point |^ as well as by a low temperature 
mean-field theory ||^, in which it was claimed that, for the four-state model {p = 4), 
only the phases (12^=), (l2'=12'=+i), (2^=3), (2'=32'=+i3), (4), and (oo) {k = 0,1,2,...) 
are stable at low temperatures, {ui . . . Ur) is a shorthand notation for the phase with a 
period consisting of r bands with ui,U2, ■ ■ ■ ,Ur layers with spins n, n + 1, . . . ,n + r (all 
modulo p) respectively. The phase (12), for example, is given by the layer sequence 



112 3 3 11 
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The ferromagnetic and chiral ground states are denoted by (oo) and (1) respectively. 
McCuUough p investigated the phase diagram for p = 3, 4, and 5 using the mean-field 
transfer-matrix (MFTM) method. From the numerical extrapolation of the data it was 
concluded that the low temperature phase diagrams for p = 3 and p = 4 were consistent 
with the results of the low-temperature series expansion @, 0]. It is interesting, that, 
for p = 5, new phases not predicted by the low-temperature series expansion were 
found to be stable at low temperatures. 

Scholten and King presented Monte Carlo simulations of the four- and the six-state 
models. They investigated especially the transition from the modulated phases to the 
ferromagnetic phase (i.e. A < |). As it was not possible to resolve particular phases, 
they determined the " interface spacing" as the average number of layers in a band for 
a given phase. They claimed that, for A = 0.45, the results were not inconsistent with 
the predictions of Yeomans. In the case p = 4 and A = 0.2 new phases with an interface 
spacing larger than the interface spacings of the phases predicted in [|| were observed 
close to the transition to the ferromagnetic phase. 

Recently the four-state chiral clock model was shown |Q to be a special case of the 
Double Ising Spin (DIS) model [jlO|, [ll|, 0, which was introduced to describe uniaxially 
modulated ferroelectrics. 

In the following new results for the four-state chiral clock model are presented. In section 
and 1^ we will reexamine the low temperature phase diagram and discuss discrepancies 
with previous results. In section ^ it is shown that the transition from the modulated 
phases to the paramagnetic phase belongs to the universality class of the 3d-XY model, 
and in section 5 short conclusions are given. 



2. The low temperature series expansion 



The present series expansion technique for the four-state chiral clock (CC4) model is 



similar to the method developed by Fisher and Selke ||13[ for the Axial Next Nearest 
Neighbour Ising (ANNNI) model. At low temperatures the reduced free energy per spin 
/ 



j^f^^rp {N is the total number of spins) may be expanded in the form [0] 



f 
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n>l 



AZjj*^ is the total contribution to the partition function from configurations in which 
n spins have flipped (as compared to the ground state). Eq, the ground state energy 
per spin, can be expressed [g, |^ in terms of the structural variables Ik = L^jL (L^: 
number of fc-layer bands; L: total number of layers): 



^0 ({4}) = --q±Jo - Ji 



k>l 



with Ji = J cos (f and 6 = tan (f — 1. The number of nearest neighbours in the 
layers is q±; it is 4 for the primitive cubic lattice. 

The contributions AZ^^ are expressed in terms of the elementary Boltzmann factors 
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with Kq = jQ/{kBT) and K = J/iksT). The reduced free energy per spin can be 
expanded in a convergent power series of w, provided that x ^ w, i.e. if Jo is large 
compared to J (which is assumed throughout this paper). The weight w results from 
changing an in-layer bond between spins with equal values to a bond between spins with 
values differing by 1. The lowest orders involved are w'^^ (overturning one spin), w^'^^~^ 
(overturning two neighbouring spins in one layer) and w'^'^^ (overturning two spins not 
being in-layer nearest neighbours). 

There are three possible environments of a given spin (the numbers in parentheses are 
the values of the spins in three consecutive layers where the considered spin belongs to 
the middle layer): (a) spins with two ferromagnetic bonds in the chiral direction (e.g. 
000), (b) spins with one ferromagnetic and one chiral bond (e.g. Oil), and (c) spins with 
two chiral bonds (e.g. 012). 

Let us discuss, as an example, the contribution to AZ^"* [first-order term in equation 
(0)] for case (a). By overturning one spin, three different final states can be obtained 
(m being the initial state): (m + 1) mod 4, (m + 2) mod 4, and (m + 3) mod 4. This 
leads to the Boltzmann factor 



^exp[- {Ef{n) ~ E,) / (ksT)] 
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It is obvious from equation (3) that the process m — > (m+2) mod 4 does not contribute 
to the lowest order term in the expansion, as it has the same in-layer Boltzmann factor 
ly^^-L as the higher order process by which the values of two uncoupled spins change by 
1. In fact, this process of the order w^'^-l does not even contribute to the lowest order 
correction term, which is of the order w^''^"^ (flipping of two neighbouring spins in one 
layer [[T^ ). 



In reference p| the following contribution to AZ^"* for the case (a) is given: 
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A comparison of equations (4) and with equation (3) reveals that the treatment of 
the in-layer bonds is erroneous in reference 0. The free energy is written in reference 
as an expansion in terms of the (erroneous) Boltzmann factor u. As a consequence, 
contributions from different orders of the expansion are treated in reference @] as if they 
were of the same order. Thus, in our example, the term x^, resulting from the process 
m — > (m + 2) mod 4 and contributing to a higher order correction in the polynominal 
expansion in w [see equation (3)], contributes to the lowest order in p[ [see equation 
(4)]. This error is repeated for all considered spin configurations and for all considered 
p-state models {p > 4), thus leading to a wrong low temperature phase diagram not 
only for the CC4 model, but also for the generalised p-state chiral clock model with 
p > 4. One should emphasise that the treatment of the in-layer bonds is correct in the 
analyses of the CC3 model [0. 

With the correct contributions, the reduced free energy (2) in first order is given by 



Ki - -K16 - {l + xy)w 
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The set of structural variables Ik minimizing / for given values of 6 and T determine 
the stable phases occurring in first order (see figure 1): the {00)-, the (1)-, and the (2)- 
phase. Phases {00) and (2) are, in this order of the expansion, separated by a boundary, 
at which all phases that are degenerate at the multiphase point and that do not contain 
1-layer bands have the same free energy. Likewise, phases containing only 1- and 2-layer 
bands are still degenerate on the boundary between the (1)- and the (2)-phase. 
One could now proceed in considering processes involving two spins, then three spins 
and so on. This is very cumbersome and only feasible for processes involving few spins. 
In the next section the phases stable in general order in the series expansion will be 
determined using a transfer-matrix method. 
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3. Transfer-matrix method 



3.1. Introductory remarks 

One should first note that the Hamiltonian (|l|) is left invariant by the transformation 
A — ^ A' = 1 - A 

ni,a — ^ K a = {-rii^a + ") mod 4. (7) 

Therefore, the phase diagram of the CC^ model is invariant under a reflection in the line 
A = |. In the following we will discuss the low temperature phase diagram for the case 
A > i, i.e. we will analyse in detail the stability of the boundary line between the (1)- 
and the (2)-phase, the phase diagram for A < | being inferred by the transformation 
(7). 

In the ground state and in the low temperature expansion every phase (y) consists of a 
periodic arrangement of a sequence of n{v) layers called i/-sequences {n{v) is the period 
of the phase]. Suppose now that in a certain order of the series expansion two stable 
phases, {pi) and (z/2), are separated by a boundary at which the phases produced by 
vi- and z/2-sequences are degenerate (see flgure |^). In first order the boundary under 
consideration separates the phases (1) and (2). At higher order a new phase {v) = (i^i^'2) 
consisting of a structure with alternating ui- and z/2-sequences might be stable in the 
vicinity of the boundary. If 

— J{u) 7 — ^— 7 — \J{iyi) 7 — T"^ 7 — \J{'^2) \°) 



is negative, the new phase has a lower free energy than the phases (z/i) and {1^2) 0; 
and it will be stabilized in the vicinity of the (ui) : {1/2) boundary (see figure 2a). The 
stability of the boundaries between the phases {ui) and (z/iZ/2) and the phases (i^i^'2) 
and (1/2) must then be examined at higher orders. If, on the other hand, is positive, 
the phase {vii'2) (and therefore every phase consisting of z/i- and i/2-sequences) has a 
higher free energy than either (z^i) or (z/2). The boundary is a true phase boundary 
which remains stable in all orders of the low temperature series expansion (see figure 
2b). 

The reader is referred to references [jl3| and |^ for details concerning the construction 
of the series expansion to general order. 



3.2. Formulation in terms of transfer matrices and vectors 

The sign of a^, and therefore the stability of the phase (z/), is determined by the leading 
term in its expansion in terms of w. This term is obtained by considering all hipping 
processes involving a spin chain of n{v) — 1 spins in n{v) — 1 different layers [0. Besides 
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the linear configuration with all ?t,(z/) — 1 spins connected, the various decompositions 
of this configuration into 2, 3, . . ., n(z/) — 1 different parts must be taken into account. 
The contributions from these processes can be written as a product of transfer matrices 
and vectors. The matrices describe a bond between two flipping spins, the vectors an 
initial or a final bond preceding the first or following the last flipped spin respectively. 
Every spin can flip to three different values and hence 3x3 matrices occur. As we 
are only interested in the sign of the we can restrict ourselves to the two processes 
contributing in lowest order, thus excluding the process m — > (m + 2) mod 4 only 
relevant for the correction term. Of course, if one considers all possible processes (i.e. 
3x3 matrices), the leading term is identical to the term obtained by the 2x2 matrices. 
This has already been noticed in the low temperature analyses of a six-state clock 
model with competing axial nearest and next-nearest neighbour couplings [|14|, where 
the corresponding 2x2 matrices have been considered instead of the general 5x5 
matrices. 

As two axial next nearest neighbours are either coupled by a ferromagnetic or by a chiral 
bond, only two different matrices are to be constructed. For a ferromagnetic or a chiral 
bond between two spins in the layers a and a-l-l one obtains, respectively, the transfer 
matrices 



and 



^-^^"U(i-^) 1-1/ r ■ ^ ^ 

The matrix elements are the Boltzmann factors for a simultaneous change of the values 
of the two spins. The first (second) row corresponds to a change Arij^Q = +1(~1) cind 
the first (second) column to Ani^a+i = +1(~1)- Every element of the matrices F and C 
is a sum of two terms, the first term resulting from changing the values of two axially 
coupled spins. As already mentioned, disconnected pairs of spins (i.e. two spins that are 
not neighbour to each other but neighbour to an unchanged spin) also contribute to the 
partition sum. Since every disconnected pair must be associated with a minus sign |13], 
the corresponding Boltzmann factors enter the different matrices with a negative sign. 
The factor w'^^ resulting from changing the in-layer bonds in layer a+l is common to 
all elements of the matrices F and C. This is a direct consequence of the fact that only 
flipping processes m — >■ (m ± 1) mod 4 are to be considered for obtaining the leading 
order in the expansion of a^. For the full 3x3 matrices this is not the the flipping 

process m — > (m + 2) mod 4 has the in-layer Boltzmann factor w'^'^^. In reference 
the phase diagram has been determined to general order using 3x3 transfer matrices. 
Due to the erroneous treatment of the in- layer interactions (see section |^) the "common 



term" uj'^^ has been factorized, thus leading, again, to the treatment of terms belonging 
to different orders as being of the same order. 

A spin at the end of the spin chain is neighbour of an unchanged spin. To determine 
the contributions of these spins, four different cases are to be distinguished: (a) the 
considered spin is the first spin of the chain and its bond to the left (i.e. to an unchanged 
spin) is a ferromagnetic or a chiral bond (subscripts / and c respectively) or (b) it is 
the last spin of the chain and its bond to the right is a ferromagnetic or a chiral bond. 
The Boltzmann factors for the flipping of these single spins are written as vectors: 



1 

y 2 \ 1 

X2 \ 1 
_ 1_ 

X 2 

1 



1 



x^w^-L (in 



y^w'^^ (12) 



b^ = j (13) 

^ = ( 5' ) (14) 

The vectors (13) and (14) do not include the Boltzmann factor resulting from the change 
of the in-layer bonds. This factor has already been included in the matrix describing 
the overturning of the two last spins in the spin chain. 



3. 3. Derivation of the low temperature phase diagram 

With the matrices (|]) and ([lOD and the vectors (11)- (14) it is now possible to compute 
the leading order term by of the quantities (and, thus, to determine the sign of a^) for 
all phases degenerate at the multiphase point and containing only 1- and 2-layer bands. 
All considered phases can be viewed as periodic arrangements of spin sequences with a 
1-layer band as the first and a 2-layer band as the last band in the sequence 0. The 
sequence z/ obtained by stripping the original sequence v by its last and first band is 
called core. All sequences based on the same core D enter in the computation of the h^: 
The sequences lz/2 and 2i>l contribute negatively, the sequences lz>l and 2z>2 contribute 
positively |jl3[. The expressions hy for different families of phases are summarised in 
table 1. 

3.3.1. Stability of some series of phases For the series of phases (l2^'^ the expression 
(see table 1) 

fei2^ =-(af-aJ) (£E)'"'g(b^-b,) (15) 

gives the leading contribution to ai2k- The four different sequences based on the core 
z/ = 2^^^ yield the four different contributions to 6i2fe. 
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The eigenvalues exp (— F-i-) of the matrix C F are real and positive. Expression (p!5|) can 
be written in the form 

6i2fc = exp + A„ exp 

with r_|_ < r_. A close examination reveals that for finite temperatures < 0, A_ > 
and A^ + A^ < 0. Thus, 612'= is negative for all fc, i.e. all phases of the form ^12'^^ spring 
from the multiphase point and have a finite stability range at temperatures above zero. 
The leading order contribution for the phases (1^2) is 



bik2 = - (aj - aj) (b^ - b,) . (16) 
The eigenvalues of the matrix C are complex conjugate. They are written in the form 

ei±i6 = exp(-ro±ir]) (17) 

with To = -i In {^f + Q > 0, n = arctan (|) and ^ = 1 - x^+^, ^2 = (1 - 
We then obtain the expression 



x) X^^". 



V2 = -(Ai + iA2)(ei + i6) + (^i-i^2)(ei-i6) = -|A|exp(-A;^o)cos(A;^] + 0) 

with I A| exp icp = Ai + iA2. The temperature-dependent quantities | A|, (p, Tq, and Q do 
not depend on k. 

hik2 is negative for small values of k. If k exceeds the value kmax = ^ (f ^ 'Z')' then 
hik2 becomes positive and, thus, all phases with k > k^ax are unstable at the considered 
point of the phase diagram. Since kmax — ^ 00 for T — > 0, there is, for every k, 
a temperature below which the phase (l''2^ is stable. Thus, all phases (^1^2^ spring 
from the multiphase point, but the higher commensurate phases disappear at higher 
temperatures. Such a cut-off of the high commensurate phases at finite temperatures is 



also observed in the ANNNI model [15 



Following the general line we also examined the series of phases <|l2'^12'^'*'^y and 
^]^fc2]^fc-i2^_ For the case <^12'^12'^+^^ we find that all these phases are stable at finite 
temperatures in the vicinity of the multiphase point with no cut-off for the phases with 
a large value of k, i.e. the results for the series (^12*^12*^"'"^^ resemble the results for the 

series ^12^^. Analysing the leading contribution for the phases ^1^21^~^2^ we find a 
behaviour similar to the behaviour of the phases <^1^2^, i.e. all phases with k < k^ax 
(the value of kmax being series-dependent) are stable and kmax — > 00 as T — > 0. 



3.3.2. Phases containing general sequences of 1- and 2-layer hands In the following 
we will show that all phases consisting only of 1- and 2-layer bands and obeying 
the rules of the structure combination spring from the multiphase point, the higher 



10 



commensurate phases of some series becoming unstable at higher temperatures. The 
leading contribution to for all these phases is of the form (see table 1) 

&^ = - (af - aj) S£ (b/ - b,) (18) 

where H is a product of powers of matrices Q_ and (GE)- The contributions of the 

first and last band are given by {s^ — aj) and Q.{hf — be) respectively. A 1-layer band 
in the core contributes a matrix C, whereas a 2-layer band yields the matrix product 
(^CF). The product over all bands in the core yields the matrix D [see equation (|18|)]. 
The diagonal elements of the matrix 

\y\^-x) -xy [l-y) [l+xy)[l-y)[l-x) ) 

are positive whereas the non-diagonal elements are negative, since y = x^^^ with x <^ 1. 
We now follow reference and introduce the unitary matrix 




All elements of U (CE) H, and therefore of U (CE) H, are positive. This is also the 

case for the two vectors (^aj — aj) ^ and (hj — b^), thus [see equation (|15|) and 
table 1] 

(aj - aj) {QrY'^ Q {hf - b,) = (aj - aj) HE {QkY'^ RRQ {h - k,) 

is positive, i.e. bi2k < 0, in agreement with the aforementioned calculations. 

Phases of the series {^12'^ 12'^"*"^^ contain a single 1-layer-band in the core yielding 

the matrix product (GE) Q (£E) with positive diagonal and negative non-diagonal 
elements for small x. Hence, the product (see table 1) 

(a^-aj) (gE)^£(gE)^g(b,-b.) 

= (af - aj) {qkY'^ (£e) q (ge) u 



is positive, showing the stability of the phases <^12'^'12*''+^^. Following this line of thought 
one easily shows that all phases appearing between the phases (2) and (12) (i.e. phases 
with only isolated 1-layer-bands in the core) are stable in the vicinity of the multiphase 
point. Indeed, as no new matrix products show up in the computation of the different 
bi,, all these expressions can be written, using the matrix as a product of vectors and 
matrices having only positive elements. 

For phases containing consecutive 1-layer-bands in the core the following additional 
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vectors and matrices may contribute to the by as can be seen from table 1: 
UQ!' (b/-b,), (a^-aj) G^^, and H {QE) Q^^R wi^^ k > 2. Introducing the 
eigenvalues of the matrix Q [see equation (p!7D], we obtain 



Vl 
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exp {-kTo) 
exp {-kTo) 



r £ 

■ _S 
X 2 



1 — x) COS kQ + X 
l-x^+^) cosm 



1+^ 

2 



1 — X 

1+5 n 



1+5 



sin fcfi 
a;) sin kQ 
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for the components of the vector v = UC'^ (h_j — b^^. Whereas vi is always positive, 
V2 > only if the inequality 



tanfcf^ < 



X 2 



(l-a:l^ 



x] 



(1 

holds. This is the case for temperatures smaller than an upper limit which depends 
on 6 and k. In a similar way one shows that for temperatures smaller than some 
fc-dependent temperature all the components of the vector (^a^ — aj ) C'^U and of 

the matrix U (CF^ C^U are positive. The free energy differences a,^ for all phases 
containing consecutive 1-layer-bands in the core are therefore negative below a certain 
temperature, i.e. these phases possess a stability region below this temperature. 



3.3.3. Conclusion The results obtained so far can be summarised as follows: All 
phases consisting only of 1- and 2-layer-bands, that can be formed by means of the 
aforementioned structure combination rules, spring from the multiphase point, where 
they are degenerate. The higher commensurate phases of some series, i.e. those phases 
formed in higher orders of the combination process, disappear again at temperatures 
individually depending on the series under consideration. 

From these results the complete low temperature phase diagram of the CC4 model 
is deduced by applying the transformation (7). At non-zero temperatures all phases 
appearing between the phases (12) and (3) are stable since the transformation (7) 
transforms phase (12) into (3) and leaves the phase (2) invariant. Some of the long 
commensurate phases appearing between the phases (1) and (12) for A > | and between 
the phases {00) and (3) for A < | are unstable at a given temperature. Upon reducing 
the temperature, more and more of these phases become stable, and in the limit T — > 
all phases obeying the rules of the structure combination are stable. Therefore, the CC4 
model exhibits a complete devil's staircase in the low-temperature limit. 



3.4. Comparison with other work 

The low-temperature behaviour of the general p-state chiral clock model was analysed in 
reference using a series expansion technique similar to the one presented here. Due to 
the incorrect expansion (see section 2) only some specific families of phases were shown 
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to possess a finite stability region at small temperatures. Especially, it was claimed that 
the phases (l^T^j with k > 2 are not stable at low temperatures, implying, due to the 
transformation (7), that for A < | a direct transition from the ferromagnetic phase 
to the (4)-phase exists. In order to corroborate these calculations a low temperature 
mean-field analyses of the CCp model was presented in |^ were it was claimed that in 
the vicinity of the multiphase point the mean-field approximation yields the same stable 
phases as reference 0. In that work the model in mean- field approximation was mapped 
onto an one-dimensional array of interacting domain walls. This mapping was derived 
under the approximation that the mean-field average spin ((cos \ni^a)MF, (sin \ni^a)MF) 
in each layer (layer spin) does only deviate from the T = value in amplitude but not in 
phase. In a detailed analyses of the mean-field phase diagram of the CC^ model Siegert 
and Everts [l^ showed that this approximation leads to a wrong phase diagram at low 
temperatures, thus concluding that the layer spin must also be allowed to deviate in 
phase from its ground-state value. This should not only be the case for the three-state 
but also for the general p-state model. The results of reference [^] for the mean-field 
low temperature behaviour of the CCp model must therefore be considered with care. 
As we have shown in the preceding sections the results of the series expansion in 
reference are erroneous due to wrong Boltzmann factors for the in-layer bonds. 
In fact, the four-state model exhibits in the low-temperature limit a complete devil's 
staircase. Furthermore, it results from our calculations that no direct transition from the 
ferromagnetic to the (4)-phase exists as phases with longer periods are stable between 
these two phases. 

In the Monte Carlo simulation of the CC4 model long-period spin patterns were 
observed when going from the ferromagnetic phase to the modulated phases at rather 
high temperatures. In view of the present work one must interpret these patterns as 
reflecting the existence of phases springing from the multiphase point and intercalating 
between the ferromagnetic and the (4)-phases. 



4. The critical behaviour 

The critical behaviour of the general p-state chiral clock model at the transition to 
the paraphase is an interesting topic since for p = 2 the chiral clock model reduces to 
the anisotropic Ising model, for p = 00 it corresponds to the classical Sd-XF model. 



Siegert and Everts JT^ showed that the CC3 model belongs to the universality class of 
the 3d-Xy model. On the basis of his contradictory MFTM results, McCullough 
speculated about a change in the universality class from 3d-Ising behaviour to Sd-XF 
behaviour for p close to 5. In the following we will show that for p = 4 an effective 
Ginzburg-Landau- Wilson Hamiltonian can be derived which can be transformed to the 
effective Hamiltonian of the 3d XF— model. 
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For the case p = 4 the Hamiltonian (|1|) can be rewritten in the form 

i Q « {ij) 

where we introduced the spin vector Sj q, = {cos ^Ui^a, sin ^nj.Q,^ and the rotation matrix 
R(A) = 



cos f A sin f A 
-sinfA cosfA 



Rotating all spins in layer a by the angle |aA, i.e. introducing new vectors 
R(aA) SjQ,, leads to the expression 



Z = ^ exp 

{2} 



1 ^ ^ 

2 51 " :i P ^ln 

ij a 13 K=l 



(20) 



for the partition function, n labelling the two spin components. The elements Ki^j/s 
of the coupling matrix are zero unless the lattice sites {i,a) and are nearest 

neighbours. 



Expression (^) may be transformed [18, IS] to 

1 



Z = CT.(l[I[Jdhi] exp 



Here C is a numerical constant, is the total number of lattice sites and I is the N x N 
identity matrix. The matrix L is given by L = /U I — K where the positive number /i is 
chosen large enough to ensure that all the eigenvalues of L are positive. 
The sum over all states can be easily computed: 

Eexp(EEE^^a<J 

{a} \ i a K ) 

= 2^ n [cosh {h\^ c„ - s„) + cosh (h\^ + c„) 



(21) 



with Cq, = cos {^al^ and = sin (^|aA^. Using the expansion 

22n-l (^22n _ 1) 



In coshx = — ^ (~1) 



n(2n)! 



{Bn- Bernoulli number) the expression on the righthand side of (21) can be written as 
the exponential of a sum of powers of h}^ and /if^: 



Ci exp 



2n k 2n—k 

EEEEc(^>^,M')EEcL+'^ 

L n k=Ql=Q l'=0 i 



J.+1' 2n-l-l 



2n-k 



(22) 
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where ') is a number depending on n, k, I and 



-1 



Introducing new variables 
the wavenumber representation leads to the partition function 



h'^ a, taking the continuum limit and turning to 



exp 



with the effective Hamiltonian 



H 



BZ (27r)^ 



-u 



T 

d^qd^q' d^q" 



BZ 



(2vr) 



Z{q) Zi-q) 



z{q) z{q)] \z{q") z{-q -q - q 



(23) 



with Tq = (Tq,Tq^ and 0'' (r) = /g^ exp(iq ■ r) r^. The integration is over the 
first Brillouin zone with q = (^qj^,qj, its components q^^ and q^ being perpendicular 
and parallel to the direction of the modulation respectively, r = (1 — 2To — T) with 
T = and Tq = varies linearly with temperature. 

In deriving equation (23) we neglected fourth and higher harmonics, i.e. fast oscillating 
terms containing exp (^in|aA^ with n > 4. Furthermore we did not include terms of 



higher than fourth order in r. If we rescale gy in the effective Hamiltonian |^ we end 
with the effective Ginzburg-Landau- Wilson Hamiltonian of the 3d-Xy model. 



5. Conclusions 



A low temperature series expansion technique is suitable to obtain exact results on the 
low temperature behaviour of the four-state chiral clock model. All phases degenerate 
at the multiphase point (T = 0, A = i) and obeying the structure combination 
rules spring from the multiphase point. Some of these phases disappear at higher 
temperatures. In the low-temperature limit the CC4 model exhibits a complete devil's 
staircase. Differences in the low temperature phase diagrams derived in the present and 
in a previous publication can be traced back to an inconsistency in the series expansion 
of the latter. Long-period spin patterns derived in the present paper as stable phases 
between the ferromagnetic and the (4)-phase and not occurring in the analyses presented 
in P], were recently seen in Monte Carlo simulations just above the boundary of the 
ferromagnetic phase. 

Furthermore, the critical behaviour at the boundary between the paraphase and the 
modulated structures follows from the derivation of an effective Ginzburg-Landau- 
Wilson Hamiltonian. It is shown that the latter can be transformed to the effective 
Hamiltonian of the 3d-Xy model. The four-state model thus belongs to the universality 
class of the XY model. 
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Tables and table captions 



Table 1. The leading ordc'rs h„ in the expansion of the quantities determing the 
stability of different families of phases consisting of 1- and 2-layer-bands. 
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Figure captions 



Figure 1. Schematic phase; diagram showing, for given small value of T, the phase 
sequence on a line along which 5 varies. Lower horizontal line: sequence in zeroeth 
order of the expansion (exact for T = 0); upper line: first order results. (5 = is the 
multiphase point. Dashed lines indicate boundary lines at which an infinity of phases 
are degenerate. 



Figure 2. Schematic phase diagram in the vicinity of the boundary between two 

stable phases {i/i) and (j^2}- The horizontal lines correspond, as in figure 1, to a given 
value of T and present results of nth and mth order, (a) The phase {viv^) is stable 
at higher order leading to new boundary lines, (b) The phase {v\V-2) is not stabilized 
leading to a true phase boundary between the two phases {v\) and (1/2). A dashed 
line indicates a boundary at which an infinity of phases are degenerate, a solid line 
indicates a true phase boundary. 
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